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The recent developments of the phase string theory for doped antiferromagnets will be briefly 
reviewed. Such theory is built upon a singular phase string effect induced by the motion of holes 
in a doped antiferromagnet, which as a precise property of the t-J model dictates the novel com- 
petition between the charge and spin degrees of freedom. A global phase diagram including the 
antiferromagnetic, superconducting, lower and upper pseudogap, and high-temperature "normal" 
phases, as well as a series of anomalous physical properties of these phases will be presented as the 
' self-consistent and systematic consequences of the phase string theory. 
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. I. INTRODUCTION 

1 Since the discovery of high-T c superconductivity two decades ago, a great effort has been put into the search for 
a unified microscopic mechanism for both superconductivity as well as anomalous spin and charge properties in the 
1 "! cuprates. 

In this paper, I shall review a systematic endeavor along a particular line approaching the doped Mott insula- 
tor/doped antiferromagnet, which has been proposed 1 as the unique property of the cuprates due to the strong 
Coulomb interaction. It has gradually become a consensus that the doped Mott insulator physics holds the key to 
understanding the cuprate superconductor, and distinguishes the latter from a conventional BCS superconductor. 

The present line, known as the phase string theory, will be characteristically different from the main-stream approach 
to a doped Mott insulator. The latter has been mainly built on the resonating valence bond (RVB) pairing of fermionic 
I i spins and the spin-charge separation in the early proposals, 1-3 the latest developments of which have been summarized 
and reviewed in Rcfs. 4,5 . 

The phase string theory, on the other hand, has been built on a singular nonlocal effect hidden in a typical doped 
Mott insulator, e.g., the t-J model. This effect is largely omitted in other approaches, but is critical in constructing a 
self-consistent theory of the doped antiferromagnet evolving continuously from the half-filling Mott-antiferromagnetic 
insulator. I shall present a rich phase diagram with complex phenomena as the physical consequences of the phase 
string effect. One will see that the RVB and spin-charge separation concepts remain essential, but they will also 
acquire a distinct mathematical characterization in this theory. 

Although the leading order results of the phase string theory share some striking similarity with the experimental 
measurements in the cuprates, in this review, a direct comparison with experiment will be minimal, partly due to 
the length limitation and partly because I wish to emphasize that the physical consequences naturally flow from 
^1 , the theoretical structure, not from a phenomenology based on the experiment. I will not be able to discuss many 
interesting theoretical efforts along different lines of thought in this brief review, also due to the space limitation, 
| which may be found through our original papers in the references. 

The rest of the paper will be organized as follows. In Sec. 2, I will discuss some important general properties of 
the t-J model, including the Marshall sign rule, bosonic RVB description, phase string effect, and exact phase string 
formalism. In Sec. 3, based on the phase string formalism, I will describe how an effective theory, known as the phase 
string model, is constructed. In Sec. 4, the physical consequences of the phase string model will be given which cover 
the superconducting phase, lower and upper pseudogap phases, high-temperature normal state, as well as low-doping 
antiferromagnetic state. Finally, a synthesis and perspectives are presented in Sec. 5. 
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II. GENERAL PROPERTIES OF THE t-J MODEL 



The t-J model is defined by H t -j — H t + Hj\ 

H t = -t^24 a c j(T + h.c. (1) 

M 

and 

H J =jJ2(s i -S j -^i) (2) 
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with a restricted Hilbert space under the no double occupancy constraint 

E c *W<i. (3) 

Due to the no double occupancy constraint, the t-J model describes a Mott insulator at half filling with 
rii = ^2 a c\ a Ci a = 1, where the hopping term H t — and the superexchange term Hj depicts the Heisenberg 
antiferromagnetic (AF) interaction in the unfrozen spin degrees of freedom. Away from the half filling, charge carriers 
are introduced with removing or injecting electrons into the system, which are known as the hole or electron-doped 
Mott insulators. Since the half-filling case is an AF spin state, the doped system can be also properly called a doped 
antifcrromagnet . 

For a bipartite lattice with only the nearest-neighboring (nn) hopping and superexchange couplings [denoted by 
(ij) in Eqs. (1) and (2)], there exists a particle-hole symmetry and one may only focus on the hole-doped side without 
loss of generality. Note that the next nn hopping may be important in real materials and its effect will be commented 
on in Sec. 5. 



A. Half-Filling: A Bosonic RVB Description 

At half-filling, the t-J model reduces to the Heisenberg model, where the physical properties are fairly well under- 
stood, as a consensus, in contrast to the doped case. This model predicts a long-range Neel order in the ground state 
as well as low-lying spin-wave excitations, consistent with experiment. 

It is important to point out that in the study of the doped case, one needs to have a correct description of spin 
correlations not only at long distance, but at all ranges. This is because the hole hopping in H t involves the nn 
sites, which will be generally quite sensitive to the short-range spin correlations. Thus, as a starting point, a precise 
description of both long and short range spin correlations at half-filling is essential. 

1. Marshall sign rule 

As proven by Marshall, 6 the ground-state wave function of the Heisenberg model for a bipartite lattice must be real 
and satisfies a sign rule. This sign rule dictates that the flip of a pair of antiparallcl spins at two opposite sublattice 
sites will always be accompanied by a sign change in the wave function. 

This Marshall sign rule may be easily understood as below. Define a spin basis state with the built-in Marshall 
signs as 

\4>) = (-1)^1 T • • • IT • • • I) (4) 

where N^ A denotes the total number of down spins at the A sublattice such that the aforementioned Marshall sign 
rule is always satisfied. Then it is straightforward to verify that the matrix element of Hj under the complete set 
{\<fi}} is negative definite 

tf\Hj\4>) < (5) 

so that for the ground state 

l*o> = 5>*l*> (6) 

{0} 

the wave function \^ is always real and positive (except for a trivial global phase). 

2. Liang-Doucot- Anderson wave function 

While the exact ground state of the Heisenberg model in two dimensions (2D) is not known, the best variational 
state proposed by Liang, Doucot and Anderson is given by 7 

l*0>= E W hh---W injn (iljl)...(i n jn) , (7) 
lEAjGB 
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where (ij) stands for a singlet spin pairing at two opposite sublattice sites i and j, and Wij, the positive weight factor 
associated with it. Since each singlet bond (ij) automatically satisfies the Marshall sign rule, the total wave function 
can be easily shown to obey the sign rule as well. 

The variational wave function in Eq. (7) can produce 7 an unrivaled accurate ground-state energy (— 0.3344J per 
bond as compared to the exact numerical value of — 0.3346J per bond for the Heisenbcrg model). Since the energy of 
the Heisenbcrg model is directly related to the nn spin-spin correlation, a good variational energy also means a good 
description of short-range spin correlations. On the other hand, this wave function possesses an AF long-range order 
(AFLRO) in the spin correlation function with a similarly accurate magnetization as the system size is extrapolated 
to infinity. 7 



B. Doping: Phase String Effect 

Now let us consider doping. In the above we have seen that the ground state wave function of the Heisenberg model 
satisfies the Marshall sign. In fact, such a Marshall sign rule would hold even at arbitrary doping, if holes remain 
static on lattice sites. Consider the single hole case for example. Define the following spin basis incorporating the 
Marshall signs similar to Eq. (4) 

|^;(n)) = (-l)<|T---| T---|) (8) 
where n denotes the hole site. Then it is straightforward to check that 

(<t/;(n)\Hj\4>;(n))<0 (9) 

still holds to ensure the Marshall sign rule. 

But once the hole starts to move, the Marshall sign rule will be scrambled by the hopping of the hole. Indeed, 
based on Eq. (8) the hopping matrix is given by 



(<f>;{m)\H t \<f>-,{n))=-ta m (10) 

where a m is the site m spin index in the state \<j>; (n)), and \cj>; (m)) is different from |</>; (n)) by an exchange of the 
spin a m with the hole at site n. Since a m = ±1, the hopping matrix element is no longer sign definite. 
Based on Eq. (10), a hole hops on the lattice will pick up a product of sequential + and — signs, 8 

l[a m = (+1) x (-1) x (-1) x ... 

C 

EE (11) 

where N* is the number of j spins which are exchanged with the hole during its hopping on a given path c. Because 
of Eq. (9), the superexchange interaction cannot "repair" such a phase string effect created by the nn hole hopping. 8 

For example, the single-hole propagator G(j,i;E) = (^ \c^ a G(E)c ia \^Q^ with G(E) = l/(E — H t -j) may be 

formally expressed as 



G(j,i;E)ocJ2(-^) Nt W(c;Nt,E) (12) 

in which for each path c connecting i and j, there is a phase string factor (— 1)^ weighted by W(c; N^;E) with 

W(c;Nl;E)>0 (13) 

at E < Eq, whose proof is based on Eq. (9). 8 Here E G denotes the ground-state energy when the hole remains 
static. Similarly, the total energy E^ of the state at a momentum k can be also formally expressed in terms of the 
Wigncr-Brillouin expansion as 9 

E k =E° G -^j2 eik<ri ~ rj)M ^ (14) 
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where 



M^]T(-l)^M[ C ;7Vj] (15) 

c,ATi 



with a positive semi-definite weight functional 



M[c;N l c ]>Q. (16) 



Physically, the phase string effect represents the transverse spin mismatches created by the motion of the hole on 
the AF background, in additional to the mismatches in the S z component. 8 The irreparableness of the phase string 
means that the three components of the spin defect induced by the motion of the hole cannot be simultaneously 
repaired via the spin flips in the superexchange term. It has been concluded that either the quasiparticle weight 
vanishes 8 or the doped hole is self-localized 9 due to the phase string effect in a spin background with long-range AF 
correlation (in the half- filling limit). 

One can similarly demonstrate the irreparable phase string effect in an arbitrary multi-hole case although the 
formulation will be slightly modified with the emergence of an extra sign factor (— 1)^ in additional to (— l)^ , 
where iV™ is the number of the exchanges between the holes and c here denotes the multi-hole paths. 10 

The significant role of the phase string factor (— 1)^ in Eqs. (12) and (15) is to cause strong phase interference 
between different hole paths and thus greatly influence the charge dynamics. But since the spin degrees of freedom 
are involved in N+, the spin correlations will also get simultaneously reshaped to minimize the total kinetic and 
superexchange energy. Therefore, the irreparable phase string effect will play a role to mediate nonlocal mutual 
influence between the charge and spin degrees of freedom as a new kind of interaction which emerges from the 
strongly correlated system of doped Mott insulator/antiferromagnet. 

The factor (-1)^ = ±1 is very singular, as a fluctuation in N+ by ±1 can result in a total sign change. So the 
quantum fluctuations are normally extremely strong especially for long paths, and there is no well-controlled quantum 
many-body method to directly handle such an effect. Alternatively a unitary transformation e* e can be explicitly 
introduced 8 to precisely keep track of the phase string effect. Then, with the Marshall sign basis \(j>) [cf, Eq. (4)] being 
changed to \(f>) — e lB \4>), the new ground-state wave function x<j> m I*g) = S{</>} X<j>W) should become more or less 
"conventional" as the singular phase string effect is now sorted out into \4>). 

The t-J Hamiltonian in this new representation, known as the phase string formalism, 8 is expected to be pertur- 
batively treatable as the singular part of the phase string effect is "gauged away" by the unitary transformation. In 
the following, we shall present such an exact reformulation of the t-J Hamiltonian at arbitrary doping. 

C. Phase string formalism 

The phase string formalism is equivalent to "bosonizing" the electron operator as follows 8 

c ta = h\b la e l ^° (17) 
where "holon" h\ a and "spinon" operators are both bosonic fields, satisfying the constraint 

cr 

The nonlocal phase string factor e lBi(T in Eq. (17) is denned by 

e <6 te = e ij[*;-*j-„*f] ( a fn_ a y t (19) 

where 



l^i \ a J 



(20) 



$° = ]Tlmln (zi-zi) , (21) 
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and 

= ^Imln {zi-zi)n? , (22) 

in which n\ a and nj 1 are spinon and holon number operators respectively, at site I with zi = xi + iyi a complex 
coordinate on the lattice. 

It is easily verified that the fermionic statistics of Cj CT is automatically ensured by e* eiCT , in which the factor (<r) Nh 
(Nh is the total holon number operator) guarantees anticommutation relations between opposite spins, and the 
factor (—1)' = ±1 (for i e even/odd) is a staggered factor added for convenience. Furthermore, the equality (18) 
replaces the original no double occupancy constraint (3) imposed on the electron operator. Therefore, this is an exact 
representation of the electron operator in the Hilbert space constrained by the no double occupancy condition. 



1. Nontrivial gauge structure 



Rewriting the t — J model using this new electron decomposition (17), one gets 8 



H t = ~t E h\h 3 b] a b la e< A ^- A %) + h.c. , (23) 

and 

^ = -^E {^iif ^ii ' ( 24 ) 

(ij) 

with 

a:, E< " ,r h - h - ■ • ( 25 ) 

Here the three link fields defined on the nn sites are given by 

(26) 



4^i^Imln[|5f], (27) 

and 

4 ee \ Im In [f^K . (28) 
The strengths of these link fields can be obtained as follows 

E^^EKt-O- ( 29 ) 

and 

for an arbitrary closed loop c such that the fluxes enclosed, ^ c A\- and ^2 c A^j, are determined by the number of 

spinons and holons respectively, in the region S c enclosed by the loop c. Furthermore, the phase </>°- describes a 
constant flux with a strength 7r per plaquette: 

E< = ±7r - ( 31 ) 
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The unique feature in the above phase string formalism of the t — J model is the emergence of three link fields: 
Afj, A^j, and (f>-. Without them, there should be no nontrivial sign problem in the Hamiltonian, because h and b 
are both bosonic fields. Namely the matrix elements of H t -j would be real and negative-definite in the occupation 
number representation of h and b. Consequently, the ground state expanded in terms of these bosonic fields would 
have real and positive coefficients, which is the case at half-filling as discussed above and the one-dimensional (ID) 
case to be discussed below. 

It is easy to see that the Hamiltonian H t -j is invariant under U(l)xU(l) gauge transformations: 

hi ^ hie**, At^Atj + iVi-tpj), (32) 

and 

b ia - b ia e^ , A% - A% + (9, - 9,) . (33) 

Thus Afj and A\- are gauge fields, seen by holons and spinons respectively, as the latter carry their gauge charges 
according to (32) and (33). 

Here Afj and A^ are not independent gauge fields with their own dynamics. Rather they are directly connected to 
the matter fields as a pair of mutual topological gauge fields. The term "mutual" refers to the fact that A\- describes 
quantized 7r fluxoids attached to the spinons, coupled to the holons. Conversely, A\- describes quantized 7r fluxoids 
bound to the holons, coupled to the spinons. 

By the construction, the phase string formalism is defined in a Hilbert space where the total S z is an eigen operator. 8 
So the total numbers of | and I spinons are conserved respectively, such that the topological gauge field A\- behaves 
smoothly as defined in (29). It is also consistent with the gauge invariance under (33). Different S z states can be 
connected by the spin flip operators, defined in the phase-string representation as 

b^H, (34) 

(a factor (— l)^ has been dropped for simplicity) and S~ = (S^)^ , and Sf = ^Z a <yb\ a b i( ,. These definitions follow 
from (17). The nonlocal phase ^ in (34) will play a crucial role in restoring the spin rotational symmetry. 
Finally, the superconducting order parameter can be expressed in the phase string representation as follows 



St = 



(-iy 



a 

= e *i(*J+*5)A9. 



with the amplitude operator given by 



(-l)'e 



3 e -m-i<K 



h\h]k 3 



(35) 



(36) 



(again the factor (— l)^ is omitted). 



2. One-dimensional case 



In the ID case, one may define 



Im In (* - Zl ) = | Q . f . > ^ (37) 

such that 

^•=^.=4=0. (38) 

Thus there is no sign problem in the phase string representation of the t-J model. It implies that the Hamiltonian 
may be treated within a "mean field" approximation. 8 Namely the holons and spinons defined in the phase string 
representation of the t-J model may be regarded as the true "free" elementary excitations. 
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However, the correlation functions will be highly nontrivial because of the singular phase string effect, which is now 
precisely kept in the phase factor of the decomposition (17) with 8 



in which 



2 

l>i 



and 

eh^E a <- ( 41 ) 

Z>i,a 

Thus, to create a hole by Cj CT , according to Eq. (39), means the creation of a pair of holon and spinon excitations 
together with a nonlocal phase shift. Denoting the average hole concentration (n^) = 5, the phase string factor in 
Eq. (39) can be rewritten as 

where 

*, = £(!-*) (43) 
is the Fermi momentum (a is the lattice constant) and 



= - 2 E - 5 ) + 2 E a < ( 44 ) 

£>i l>i,ot 

with (A$j) = 0. While the leading term of the phase string factor reproduces the correct Fermi momentum kf for 
the electron system, the fluctuations in A$j will be responsible for reproducing 8 the correct Luttinger liquid behavior 
known from the large-?/ Hubbard model. 

The important connection between the phase string effect and the Luttinger liquid in ID has been first established 
previously in a path-integral study 11 of the large- J7 Hubbard model. 



3. Two-dimensional case 



At half-filling, H t has no contribution due to the no double occupancy constraint, and under a proper gauge choice 
one may set = in Hj. In this limit, there is no nontrivial sign problem in the 2D Hamiltonian which is 
fully bosonized. This is the case previously discussed in Sec. 2.1, where a precise bosonic RVB description of spin 
correlations in all ranges of length scale is available, which can serve a very good starting point for the doped case in 
2D. 

In contrast to the full bosonization case at half-filling, as well as in the ID case, the nontrivial phase structure 
emerges at finite doping in 2D, which are represented by the link fields, Afp 4>ij, and A\-. These link phases can no 
longer be "gauged away" here and they completely capture the essential sign problem (i.e., the phase string effect) 
of the doped t-J model. These gauge fields are generally well controlled in the regimes of our interest: 0° is a 
non-dynamic phase describing a constant tt flux per plaquette; Afj is cancelled when spinons are RVB paired at 
low-temperature phases; remains weak at small doping or well behaves if the holons are coherent. Therefore, 
these gauge fields will be well tractable at least in low doping and low temperature regimes. 

It is noted that the bosonization decomposition (17) was actually first obtained 12 based on optimizing a slave- 
boson mean-field state using a "flux binding" scheme. Similar procedure has been also employed recently to get 
essentially the same bosonization decomposition in Ref. 13 . This bosonization decomposition may be also regarded as 
the mutual- semion scheme as described in Ref. 8 without explicitly breaking the time-reversal symmetry. It is thus 
clearly distinguished 12 from an earlier flux-binding construction leading to a slave-semion type of formulation, 14 or a 
variant of it in a more complicated semionic representation proposed 15 in literature. 
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D. Wave function structure 



In the above, we have discussed how the intrinsic phase structure of the t-J model can be revealed in the exact 
phase string formalism. In the following we further examine the corresponding wave function structure. 

A wave function ip e in the electron c-operator representation can be related to ip b in the full bosonic h and b 
representation of the phase string formalism by 16 

^ e (*i>" ' -!*M;ji,- • -Jnc-m) =K. i> b (ii,- ■ -,i M ;ji,- ■ -,3N e -M\h,- • -JnJ (45) 

where the f spin electron sites, {i u } = i\, ■ ■ •, Im, and the J. spin sites, {jd} = ji, • • -,jN e -M, and {1^} = l\, ■ ■ •, l N h 
denote the empty sites that are not independent from {i u } and {jd} under the no double occupancy constraint. Here 
and below, we use i to specify an ] spin, j a J. spin, and I, a holon, where the subscripts u, d, and h label the sequences 
of the | spins, J. spins, and holons, respectively. 
According to Eq. (17), the JC factor is given by 16 



JC = JG, (46) 

where 

j = n - z * u < ) n ( z id - z h ) i[( z L - z jj n k ~ ^ i n ~ z ^\ n ^ ~ z ^ \ ^ 

u<u d<d' h<h' uh dh 



and 



in which the coefficient C is given by 



G = C-\-l) N * H , Zlu _ Zlh , , (48) 

uh ^ Zlh ^ 



C=\J\= H\z k -z m \ , (49) 

with k and m running through all lattice sites such that C is a constant. 

It is easily seen that the Jastrow-like factors in J automatically enforce the single occupancy constraint: J vanishes 
if two spinons (or holons) occupy the same site, or if a holon and a spinon occupy the same site. The factor J further 
explicitly captures the fermionic statistics of the electrons. Therefore, the no double occupancy constraint, which 
has been considered as one of the most important effects but difficult to tackle with in the t-J model, is no longer 
important in the phase string representation tp b , since J in (46) naturally plays the role of a projection operator. 
This may be understood in the following way. In the phase string representation, the effect of K. in the original 
V> e , is transformed into the topological gauge fields, Afj and Afj, in the Hamiltonians, (23) and (24), which describe 
spinons and holons as mutual vortices, as perceived by each other. This clearly implies a mutual repulsion between 
two species, since a spinon cannot stay at the center of its vortex (where a holon is located), and vice versa. Thus 
the constraint that a holon and a spinon cannot occupy the same site is now reflected in the interactions present in 
the new Hamiltonian, and the condition (18) is not needed as an extra condition to enforce. Note that the constraint 
(18) also requires the hard core conditions among the holons or spinons themselves. But since both holon and spinon 
fields are bosonic fields, local hard core exclusions usually do not involve the sign change of the wave function. Hence, 
in the phase string representation, the local constraint (18) is neither crucial nor singular, as far as low energy physics 
is concerned. 

Finally, the singular phase string effect is captured by the factor Q in JC. Firstly the sign factor (—1)^ can be 
identified with the Marshall sign, and N\ denotes the total number of | spins in sublattice A (note that it is equivalent 

to the previous definition using (—1)^, by a trivial global sign factor). Then the phase factor Yl h , '"_ lh , will 

\ Z iu Z l h \ 

describe the phase string effect - disordered Marshall sign. Note that it is asymmetric with regard to t and { spins: 
it only involves an j spin complex coordinate z* and a holon coordinate z ; * ft , and then will acquire the following 

additional phase as a hole moves through a closed path c, Q — > Q x (— 1) W <= T , with the displaced spins being restored 
to the original configuration by the exchange term Hj. 16 



9 



III. PHASE STRING MODEL: EFFECTIVE THEORY 

The exact phase-string formalism of the t-J Hamiltonian provides a new starting point to construct an effective 
theory which can be smoothly connected to the better- understood half-filling limit. The gauge structure in the phase- 
string formalism is a very useful guide for such a construction as it generally cannot be spontaneously broken according 
to the Elitzur's theorem. 

A. Phase string model 

Based on the t-J Hamiltonian in the phase string formalism, a minimal effective model may be written down as 
follows 16 - 17 



-Hstring = H h + H s (50) 



with 

H h = - 



t h Y J {^ Ah+leA ^)h\h J +h.c. (51) 

(ij) 

JsT, (e^)blb]_ a + h.c (52) 



This model remains invariant under the gauge transformations, (32) and (33), and is thus a gauge model, known as 
the phase string model. 

The U(l)xU(l) gauge invariance here corresponds to the charge and spin S z conservations of the holons and 
spinons, respectively, which ensures the correct quantum numbers in such a spin-charge separation description. This 
is in contrast to the slave-boson U(l) gauge theory 18 where both holon and spinon carry partial charges. In Eq. (51), 
an external electromagnetic gauge potential Afj is explicitly introduced which couples to the holon field carrying an 
electric charge +e. By contrast, the spinon field does not carry the electric charge and thus describes a charge neutral 
and spin-1/2 object, which can directly couple to the external magnetic field B only by a Zeeman term 

H% M = -» B BJ2™1. (53) 

Note that, without loss of generality, the magnetic field will be always added along the spin quantization S z axis due 
to the requirement of the S z conservation in the phase string formulation. 
The global conditions of 

= NS, (54) 
N(l-5), (55) 



E 



n h 



can be added to Hh and H s by the Lagrangian multipliers, A^ and A, respectively. Due to the relaxation of the 
local no double occupancy constraint, to avoid the short-distance uncertainty at each center of a 7r-flux tube, on the 
right-hand sides of Eqs. (29) and (30), the distribution of a holon or spinon at site I should be understood as being 
slightly smeared within a small area centered at I. 

Based on the spin operators defined in the phase string representation like Eq. (34), It is straightforward to verify 
the spin rotational invariance of the phase string model 

[String, S] - (56) 

where S = J2 i Si, by noting that $^ — = 2A^ (using the fact that the core of each flux-tube being slightly smeared 
within a small area as mentioned above). The time- reversal symmetry at A e - = can be also easily verified by noting 

that b\ a — ► ab\_ a , h\ — > h\, A\- — > A^, and A\- — ► —A\-, according to their definitions, under the time-reversal 
transformation. 
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B. Topological gauge structure and mutual Chern-Simons description 

The phase string model is uniquely featured by the two topological gauge fields, Afj and A\y According to Eqs. 
(29) and (30), the holons in (51) feel the presence of the spinons as quantized ir fluxoids through A 8 -, which reflects 
the nonlocal frustrations of the spin background on the kinetic energy of the charge degrees of freedom. Vice versa 
the spinons also perceive the doped holes as tt flux quanta through A\- , which represents the dynamic frustrations of 
the doped holes on the spin degrees of freedom. 

It is instructive to reformulate the above phase string model in the following path-integral formalism 19 



Z = J DhDbi i DbiDA s DA h exp ^— dr J 



d 2 rL stling (57) 



in which the Euclidean Lagrangian of the phase string model is given by 

^string = L h + L s + L C s (58) 

where 

Lh = J2 h \ & - iA oW] h-th^, (e iA "h\hj + c.c.) (59) 
/ </./> 

L s = E 5 t & ™ Ah M h ™ - J * E (f° A HA-° + c - c ) ( 6 °) 

i<J (ij)o- 

Lcs = ^e^A;(I)d„A h x (i) (61) 

For simplicity, two chemical potential terms enforcing the global constraints (54) and (55) are not included in L str ing- 
In such a Lagrangian formalism, two matter fields, bosonic spinons and holons, are minimally coupled to the 
U(l) gauge fields, A 8 and A h , whose gauge structure is decided by the mutual-Chern-Simons term Lcs m (61), in 
replacement of the original topological constraints (29) and (30). So the phase string theory is also known as the 
mutual Chern-Simons theory. The time-reversal, parity, and spin rotational symmetries can be explicitly shown to be 
preserved. 19 

Note that the original constraints (29) and (30) can be obtained by the equations of motion for the temporal 
components Aft and A^: 



dL 



dA 8 (I) 
dL 



= => e al3 A a A h p(i) = ivn 1 } (62) 
= 0=J>e a "A a ^(J)= Truant (63) 



with A a Ap(i) = Afy(i + a) — Afy(i) and A a A|(/) = A s p(I) — Afy(I — a). Here the indices a and are used to denote the 
spatial components (a, (3 = x, y), and the lattice gauge fields Aj 7 = A s a (I) (J = I — a) and A^- = A^(j) (i = j + &). 
The lattice sites / and i refer to two sets of "dual" lattices where the holons and spinons live on, respectively, which 
is a technical way to "regulate" the theory at the short-distance, which is not expected to change the low-energy 
physics. 19 



C. Bosonic RVB order parameter 

To justify the above phase string model, let us first consider the superexchange term Hj in Eq. (24). 
Hj is expressed in terms of the RVB operator A*- which is invariant under the gauge transformation (33). It is 
natural to define the bosonic RVB order parameter 17 

A s = (AfA (64) 

for nn sites. At half filling, A s reduces to the well-known Schwinger-boson mean-field order parameter 20 A SB = 
b ir7 bj- a ) = ( A^ B ) as A\, = 0. Since A s ^ up to a temperature <~ J/&b at half-filling, A s defined in Eq. 
(64) is expected to survive and persist into a finite doping and lower temperature regime. 
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It is interesting to distinguish the Schwinger-boson order parameter and the bosonic RVB order parameter at finite 
doping. By using the aforementioned unitary transformation 8 A^ B — ► e* e A^ B e _ie , the Schwinger-boson order 
parameter can be expressed in the phase string formalism as 

A| B = e* (1/2) [d * {l)+0 i Afj (65) 

with 9i(l) = Im In [zi — zi). So Afj 8 and Afj differ by a phase factor which is composed of 2tt vortices with the 
cores located at the hole sites, I's, with nf = 1. Namely, each doped hole will induce a 2ir vortex in the Schwinger- 
boson order parameter. The general topological vortex excitation in the Schwinger-boson mean-field state has been 
previously discussed by Ng. 21 The case that a doped hole is bound to the vortex core has been discussed as one of 
several possibilities there. It turns out that since the bare hole will induce a highly frustrated phase string effect, its 
binding with a topological vortex will be stabilized, as the composite object can effectively erase the singular effect 
and is thus in favor of the kinetic energy. Such a composite object is nothing but the bosonic holon in the present 
phase string formalism. Generally the bosonic RVB state with A s ^ means that the original Schwinger-boson order 
parameter is phase disordered with A SB = according to Eq. (65), unless the holons are localized which corresponds 
to a low-doping AF state to be discussed in Sec. 4.5. 

Based on Eq. (64), a "mean-field" version of the superexchange Hamiltonian may be rewritten as 17 



(ij)a i \ a / 



(66) 



where 



J s = Jeff A s /2 (67) 



where J c ff ~ J is a renormalized superexchange coupling to be further specified below. The Lagrangian multiplier A 
is introduced to implement the condition (55). 

In the limit where the gauge field A^j can be treated as a static field, e.g., in the holon condensed regime, Eq. (66) 
can be straightforwardly diagonalized by a Bogoliubov transformation 17 

b t a = ( U mlma ~ v mlL-a) i ( 68 ) 

m 

, v 1/2 , J/2 

as H s = Y, m a E mllnal,na + const., where u m = -j. + lj , Vm = -i_ - lj sgn(£ m ), and E m = 



A - C m . The La grangian multiplier A is determined by enforcing (b\ a bicrj = (1 — S)N. The wave function 

w ma and the spectrum £ m are determined by the following eigen equation, 

UWmaii) = -Js e " AJiffl ™(i) • ( 6Q ) 
j=nn(i) 

Self-consistently, an another gauge-invariant mean-field order parameter 

({e^)blb 3a ) nn ^Q (70) 

in this mean-field scheme. 

The corresponding ground state is given by 16 

|RVB) MF = cxp ^ W iMAl j 1°) ' ( 71 ) 

where the RVB amplitude Wij = —J2m ^f Lw m<y{^) w rna{i)- It can be further shown 16 that ^ 0, only if i,j e 
different sublattices, i.e., the RVB amplitude only connects | and J. spins on opposite sublattices. The calculated 

\Wij\ is plotted in Fig. 1 for i and j belong to different sublattices, which shows that while oc |rij| -3 > 0, for 
|rjj| 3> a at half-filling, its absolute value becomes exponentially decay at finite doping: 16 

|p- -i 2 

\Wij\ oce—&- (72) 
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FIG. 1: The bosonic RVB amplitudes, \Wij\, at three different hole concentrations. The inset shows the result for half filling. 



with a characteristic length scale which is to be identified as the spin correlation length later 

£ = aJ = V2a c , (73) 
V ttO 

where a c = a/y/nS is the typical "cyclotron" radius decided by Eq. (69). 

At half-filling, the ground state |RVB)mf qualitatively describes the AF correlations fairly well. 20 In particular, one 
may define 

|RVB) =A|RVB) MF (74) 

where P s denotes the projection to a single occupancy with each lattice site occupied by one spin. |RVB) is equivalent 
to Eq. (7) with Wij shown in Fig. 1. Based on Eq. (74), a numerical calculation 16 gives rise to (Sj • Sj) = 
—0.3344(2) and a staggered magnetization m = 0.296(2) with the maximum sample size 64 x 64. These results are 
essentially the same as the best variational result 7 and also compare extremely well with exact numerical results of 
(Si • S 7 ->„„ = -0.3346 and m = 0.31. 

\ J > tin 



D. Ground-state wave function 

At finite doping, it is easy to see that a simple ansatz for the full ground state: \^b) ~ h\ h\ ...\W/B) does not work 
here because Eq. (70) would lead to vanishing (H t ). 

According to H t in Eq. (23), the hopping of a holon should be always accompanied by a spinon "backflow" or vice 
versa. But this is not properly "registered" in utilizing Eq. (70). It implies that the correct way to add holes in the 

ground state should be 16 \^g) ~ ••• {h\bi ia ^j {h\ 2 h 2(y ^j ...|RVB) or more precisely 

l*G> = E fhiilh}) ® n({/ h }) |RVB), (75) 

where ip h describes the bosonic holon wave function, while the "spinon backflow" operator is given by 

H ({*/,})= E Z({l h },{a h })b hai b h <r 2 - (76) 

where the summation over the spin index at is under the constraint 1/2 a h = S z (note that S z = in |RVB)). 
Here |RVB) is understood as describing the "half- filling" spin background even at finite doping. Obviously the no 
double occupancy constraint is automatically satisfied by \^>g) defined in Eq. (75) at arbitrary doping. 
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Then, in accordance with the gauge invariance under (32) and (33), in minimizing G\Ht\^ g) ■> fh wn l be de- 
termined as the ground state of the effective hopping Hamiltonian Hh in Eq. (51), while the rcnormalizcd hopping 
integral th is given by 

t h = (RVB|nt ^^e-< - +ff ^)n|RVB>. (77) 

(7 

In the dilute hole limit where the correlations between the "backflow spinons" are negligible, Z ({lh} , {&h}) may be 
reduced to a product of the single-spinon wave function Z a (l), which can be variationally determined by optimizing 
th as the maximal eigen value of the following eigen equation 16 

{-t h )Z a (i) = -\ J2 e-<-^Z a (j) (78) 

j=nn(i) 

with t = + ^ i, n b = 1 — S. Numerically th thus determined is weakly doping dependent with \th\ — 0.68t ~ 
0(t). w 

Self-consistently let us go back to check the superexchange energy since |RVB) is now modified by n|RVB) for the 
spin degrees of freedom: 

(*g\Hj\* g ) = ~E E l^l 2 (RVB|nt(Af J ) t A^n|RVB) 

(ij) {lh}¥=i,j 

- - J ~f E< RVB I ( A y) f A h IRVB) (79) 

(ij) 

Note that E {lh} ^ l^l 2 = 1 - 25 + 0(5 2 ). And (RVB|nt (A?.) f A|. n|RVB> ~ /(<5)(RVB| (A?.) f |RVB) with 
/(<5) = 1 - 2( 5 - 1)6 + 0{5 2 ), (g > 1), such that 

J off = J(1- 2g6 + 0(S 2 )) (80) 

at small doping. Here g has been empirically determined 22 by g = 2 by comparing with the experimental measure- 
ments. The mean- field treatment of the last line in Eq. (79) leads to the effective spinon Hamiltonian (52). 
Generally speaking, in order to minimize (^g\Hj\^>g): the "backflow spinons" in n are better paired up 16 



n({Z/»}) oc exp 



u'e{i h },<r 



|RVB) 



(81) 



with G a u , = Z a (l)g(l — l')Z_ a {V), where g(l — I') denotes the pairing amplitude between the two "backflow spinons", 
which is no longer restricted to the pairing between two opposite sublattices. The paring g(l — l') is expected to reduce 
the hopping integral th and enhance J e s, but the detailed values of them will not affect the general properties of the 
phase string model discussed below. 



IV. PHYSICAL CONSEQUENCES 

The minimal phase string model is composed of Eqs. (51) and (52). In the following we shall sec that such a simple 
model will possess a rich phase diagram unifying the AF state, the SC phase, the pseudogap regime including both 
the upper and lower pseudogap phases, as well as a high-temperature "normal state". The richness of this model 
can be attributed to the unconventional competition between the charge and spin degrees of freedom via the mutual 
Chern-Simons gauge structure. 
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A. Superconducting (SC) phase 

The SC state is a simple self-consistent solution of the phase string model at finite doping. 17 First, the bosonic 
holons will experience a Bose condensation at T = if A\- = A\- = in H h . Once the holons are condensed, the 
gauge field will reduce to a non-dynamic Afj to describe a uniform flux of strength 

E D 4 = (82) 

per plaquette. Then, according to H s , a gap will open up in the spinon spectrum, such that the fluctuations of A*j 
get gapped, which in return self-consistently ensures the holon condensation in Hh- 

With the holon condensation 7^ 0, the amplitude of the SC order parameter, Eq. (36), becomes finite: 

= (83) 

and in the ground state, the phase coherence 

'e-^Ml/O (84) 



can be realized because of a finite-range RVB pairing of spinons with a finite excitation energy gap (cf. Sec. 4.1.6.). 
Then the superconducting order parameter defined in Eq. (35) gains a finite mean value 

Af) + 0. (85) 
Note that the phase factor e - ^ 1 / 2 ) ( $i + * j ) will also decide the d-wave symmetry of Afp ( A° is s- wave-like in general). 23 

1. Ground- state wave function 
The holon condensation as the solution of Eq. (51) at A e = may be approximately treated as an ideal one with 

ip h (h,l 2 ,—,lN h ) = const. (86) 
Then the corresponding ground state of Eqs. (75) and (81) is simplified to 16 

r .-\N h /2 

|*g)sc = const. D | RVB) 

= P Nh cxp It)] |RVB) (87) 

in which P^ h denotes a projection onto a AT^-hole state and 

D = J2 G ii( h l b ™){ h fo-°)- ( 88 ) 



Equation (87) implies (D) ^ 0. Then, in terms of Eqs. (36) and (35) 16 



ij 

= ^E^< e - l(1/2) ^)A| c )^o 



where G t] = g(i-j) Z <J (i)Z_ a (j)(-iye l ^ +l ' l '^ +laA ^ is a s-wave constant 16 based on Eq. (78). Thus, it confirms 
again that the SC off-diagonal-long-range-order (ODLRO) [Eq. (85)] is established once the phase coherence (84) is 
realized in the ground state. 
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2. Generalized Ginzburg- Landau (GL) equation 



In the SC state, one may treat the charge condensate in terms of a slowly varying continuous field, (hi) — > ijj h (ri). 
Note that in the continuum limit Eq. (51) can be rewritten as 

H h = — !— / d 2 r hUr) (-iV - A s - eA e f h(r) (90) 
2m/, J 

where mt — (2ti l a 2 )~ 1 and A s is the continuum version of Afj — ■ A s given by 

A*(r) = \ J dV K(r') - n\(r')] (91) 

with n b a (Ti) = n\Ja 2 . 

By noting that the holons here are hard-core bosons with a repulsive short-range interaction, one may generally 
write down the corresponding GL free energy Fh = J d 2 r fh where 24 

fh = f h + a Kl 2 + ||V-J 4 + ("iV - A s - eA e ) 2 ^ h (92) 

with f? denoting the "normal state" free energy density. And the "supercurrent" density is given by 

J(r) =-^~ [^(r)V^(r) - V^(r)^Jr)] - - + — ^(r)V» fc (r) . (93) 

These equations are similar to an ordinary GL theory describing a charge +e Bose condensate coupled to an external 
electromagnetic field A e , except that ip h is further coupled to the spin degrees of freedom through the vector potential 
A s . It means that each isolated spin (spinon) will register as a ±7r flux tube in Eq. (92) to exert frustration effect 
on the charge condensate. Thus, such a generalized GL must be coupled to the spinon Hamiltonian H s to govern the 
basic physics in the SC state. 




FIG. 2: Basic energy scales, E^ pper and E g , determined by the dynamic spin susceptibility function \ (Qj w ) at T = 0. The 
peak positions of x" at 5 = is shown in the energy and momentum (along the Q x = Q y axis) by the dotted curve, which 
tracks the spin wave dispersion with Eg Ppcl ~ 2.3 J and E g = 0. The upper-bound energy _E" ppcr monotonically decreases with 
increasing doping from 0.05, 0.125, to 0.2. E g denotes the resonancelike peak energy at Qaf = (vr, n), which emerges in the SC 
state. Note that the finite horizontal bars at finite doping indicate the momentum widths for these non-propagating modes. 25 



3. Non-BCS-type elementary excitation: S = 1 spin excitation 



As outlined in Sec. 3.3, the spinon Hamiltonian H s can be diagonalized under the condition (82). Figure 2 shows 
the dispersive behavior of the 5 = 1 spin excitation based on the peak position of the spin dynamic susceptibility 
x"(Q, w) at different doping concentrations, which clearly depicts how the spin excitation evolves from the spin-wave 
picture at half-filling (dotted curve) to the non-propagating modes (solid bars) in the SC state. 25 
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FIG. 3: Left panel: x"(Qaf,w) shows a resonance peak at energy E g [Qaf = (7T, 7r)] at <5 = 0.125. Inset: the evolution of i? 9 
as a function of 8. Right panel: x"(Qap,^) with incorporating the fluctuational effect induced by the charge inhomogeneity. 
Inset: the local susceptibility Xl( w ) m tne same situation. [From Ref. 25 ] 




FIG. 4: (a) Momentum distribution of x"(Q, E g ), scanned along the diagonal direction Q = (q, q) at various hole concentrations. 
The intensities are normalized at the maximums. Inset shows that the data in the main panel can be well fit into a Gaussian 
function cxp( — (Q — Qaf) 2 /2ct 2 ), with a = \/2a/£ = \Tk8~, as shown in (b). [From Ref. 25 ] 



A clear spin gap is opened up at E g in the SC state with the gapless spin wave replaced by a resonancelike mode 
(Fig. 3) around the AF wave vector Qaf = ( 7T , 7T ), whose doping dependence is also shown in the inset of the left 
panel in Fig. 3. Such a resonance mode has a finite width in momentum which implies a finite spin correlation 
length £ [Eq. (73)] as shown in Fig. 4. Furthermore, the spatial charge inhomogeneity can affect the width of the 
resonance peak via A h as shown 25 in the right panel of Fig. 3, in which the inset illustrate the local spin susceptibility 
Xl(^) = (1/A0X)q x"(Qi w )- Finally it is noted that the envelop of the high energy 5=1 excitation still roughly 
tracks the spin wave with a softened upper bound energy £?" ppcr , which decreases monotonically with doping [Fig. 2]. 



4- Non-BCS-type topological excitation: Spinon vortex 

In the above we have examined the 5 = 1 spin excitation which is composed of a pair of 5 = 1/2 spinons according 
to Eq. (34). However, a single 5 = 1/2 spinon excitation will be "confined" in the SC state, i.e., will not appear 
in the finite energy spectrum. In this sense, the above 5 = 1 excitations are true elementary ones, which do not 
fractionalize. We shall elaborate this as follows. 

It is convenient to rewrite the nn SC order parameter as the mean value of (35) in the continuum version (without 
considering the d-wave symmetry of the relative coordinate for simplicity): 24 

A sc =A (e^W) (94) 

where the amplitude 



A = A s (VD 2 



(95) 
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FIG. 5: (a) An isolated spinon is always associated with a 2n vortex in the phase of the SC order parameter, (b) A pair of 
vortex-antivortex with spinons located at the cores. 

and the phase 

$ s (r)= J d 2 r' Im In [z - z'\ [n\(r') - nj(r')] . (96) 

From Eq. (96), it is clear that <& s describes 2ir phase vortices whose cores are centered at spinons: i.e., <£> s — ► $ s ±27r 
or A sc — > A sc e ±27 ™ each time as the coordinate r continuously winds around a spinon once according to Eq. (96). 
In other words, a spinon is always associated with a 2tt vortex in A sc , known as a spinon-vortex composite 24 ' 26 which 
is schematically illustrated in Fig. 5(a). A spinon-vortex and -antivortex pair of a finite range of separation will result 
in the cancellation of the phase $ s at a large length scale as shown in Fig. 5(b). 

By writing 

^ = sfp~ h ^ h (97) 
a London equation for the supercurrent based on Eq. (93) is given by 24 

J(r) = fh. \y<f, -A s - eA e ] (98) 
m h 

Since each unpaired spinon will contribute to § c dr ■ A s — ±ir in terms Eq. (91) if the loop c encloses such a spinon, 
a minimal supercurrent vortex centered around it is then given by 



{dr- J(r) = ±n-^- (99) 
Jc m h 



at A e = according to Eq. (98). 

1 z x r 

For a single spinon vortex centered at the origin, we have A s (r) = - — , for distances r » a c ~ £, the size of 

2 r z 

the vortex core. Using J = — j^-A s , one can estimate the energy cost of a spinon-induced vortex current based on 
Eq. (90) by 24 



E v = - J d 2 r A s • J - J d 2 r Ph 



(A 



s\2 



2m h 



Ph 
2m h 



J d 2 r(A s ) 2 (100) 



= H£h. [ dr- oc In— , (101) 

Am h J r a c 

where L is the size of the sample. Thus one concludes that a single 5 = 1/2 spinon excitation is forbidden owing to 
a logarithmically diverging energy. 24 

5. Topological defects: Flux quantization and Zn impurity 

The phase string model predicts that an isolated spinon excitation is a topological vortex which cannot live alone 
in the bulk of the SC state. So there is no electron fractionalization at low energy and long distance. However, in the 
presence of two kinds of special defects in the superconductor, a single spinon excitation can be naturally induced as 
a unique prediction of the model. 
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FIG. 6: Flux quantization occurs at hc/2e, with a bosonic S = 1/2 spinon trapped inside the core. 



a. Flux quantization In the presence of magnetic field, using Eq. (98) we get 

— Ij{r)-dr = 2nn- (f dr ■ (A s + eA e ) , (102) 

Ph Jc J a 

where the integral is over a closed loop c and n, an integer. Now suppose that the integration is carried over a loop 
that is far away from the magnetic vortex core, where J = 0. Then one arrives at 



lim-ejdv- A e j - j dr ■ A s = . (103) 

If A s = 0, the magnetic flux will be quantized at 27m in units of hc/e; i.e., the minimal flux quantum in this case is 
hc/e = $07 as expected for a charge e Bose system. However, the presence of A s changes the quantization condition 
radically. Suppose there is one excited spinon trapped in the core of a magnetic fluxoid [Fig. 6(a)]. Then, from Eq. 
(103), we obtain the minimal flux quantization condition 24 

dr ■ A e = ±tt , (104) 

which is precisely the quantization condition at </> = $0/2 = hc/2e. As the holons do not distinguish between internal 
(fictitious) and external (magnetic) flux in (92), they still perceive a total flux quantized at <I>o [see Fig. 6(a)], even 
though the true magnetic flux quantum is <fi . 

Therefore, an excited spinon or free 5=1/2 moment will be "nucleated" inside the core of a magnetic vortex, 
which plays a crucial role to ensure the flux quantization at 4>q. 24l But normally the Kondo screening effect due to 
the coupling to the background quasiparticles may complicate the analysis of possible experimental behavior of such 
a free moment at very low temperatures. 

b. Zn impurity In the t-J model, a zinc impurity is treated as an empty site with excluding the occupation of 
any electrons. If one reformulates the model, with the defect, in the phase string representation outlined in Sees. 2 
and 3, it is found 27 that such a "vacancy" site will behave like a topological defect which induces a vortex current in 
the resulting phase string model, as shown in Fig. 7(a). A heuristic understanding is to imagine exciting a spinon in 
the pure system at a given site and then freezing its spin such that the superexchange coupling with the surrounding 
spins is effectively cutoff. Nor a holon can hop to this site such that the effect of a zinc impurity is effectively created, 
which is nothing but a vortex in Fig. 7(a). 

Now it is natural to see why a zinc impurity will generally induce a spin-1/2 around it in the SC state: Such a 
Zn- vortex would cost a logarithmically divergent energy and thus must be "screened" by nucleating a neutral 5=1/2 
spinon which carries an antivortex and is bound to the latter, as shown in Fig. 7(b). 
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FIG. 7: (a) A vacancy (zinc impurity) always induces a vortex-like supercurrent response in the SC phase due to the phase 
string effect, (b) To compensate such a vortex effect, a spinon, which carries an antivortex, has to be trapped around the zinc 
impurity, giving rise to a local 5=1/2 moment. 




FIG. 8: Left panel: Spinon density distribution (nt) = ^L^so-y around the zinc impurity at 8 = 0.125. Right panel: The 

distribution of [Sf ) near the zinc impurity, with the scan along the dashed direction shown in the inset, where the zinc site is 
marked by the filled circle. [From Ref. 27 ] 



Once the topological origin of the S — 1/2 moment around a zinc impurity is established, a simple effective 
description of the system with one zinc impurity can be developed 27 based on a "sudden approximation" using the 
phase string model. The physical consequences, as clearly illustrated in Figs. 8, 9, and 10, are that the S = 1/2 free 
moment induced by the zinc impurity distributes in an AF staggered pattern around the "vacancy" site. 



6. Superconducting transition 



We have seen that a single excited spinon is a topological vortex. The interaction between spinon-vortices can be 
obtained by substituting Eq. (97) into Eq. (92) and carrying out the area integration 26 



F h = J J<fv l0 l\ 2 [E Q «<( r i)] V(r 12 )[£ p P4(r 2 ) 
in which a, /3 = ± refer to the signs of vorticities carried by spinons and 

V(r 12 ) 



+ const. 



nPh ln l r i ~ r 2| 

inih r c 



(105) 



(106) 



with r c ~ a. Eq. (105) is similar to the XY model, except that A s introduces tt instead 27r vortices and the vortex 
cores are attached to spinons which have their own dynamics as governed by H s with an intrinsic quantum length 
scale a r . 
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FIG. 9: Left panel: Distribution of the contributions to l/ 63 TiT from individual sites near the zinc impurity, at temperature 
T = 0.0067J. Right panel: l/ 63 TiT vs. T at different sites. Solid curve with crosses: from the site far from the zinc impurity; 
Dashed curve with triangles: the nn site near the zinc; Dashed curve with circles: the next nearest neighbor (nnn) site near 
the zinc; Solid curve with squares: average over all sites in a 16 x 16 lattice with one zinc. [From Ref. 27 ] 
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FIG. 10: Uniform spin susceptibility in the pure system is shown by the solid curve with full circles; The case with one zinc is 
illustrated by the solid curve with triangles; The dashed curve is a fit by x = 0.2390/T. [From Ref. 27 ] 



These spinon- vortex composites form vortex-antivortex bound pairs at low temperature [cf. Fig. 5(b)]. Such a 
binding directly results in a phase coherence 

,i* a (r) 



e»*W^0 (107) 
of the SC order parameter (94) . A detailed renormalization group analysis leads to the following T c relation 28 

T ^W, (108) 



which connects the phase coherence temperature with the spin resonancelike energy E g (cf. Fig. 3). 

The superconducting phase coherence (107) implies that spinons are confined in the bulk where a single spinon- 
vortex excitation costs a logarithmically divergent energy. In this case, a finite energy elementary excitation is an S = 1 
spin excitation composed of pairs of spinon-vortices, whose minimal excitation corresponds to the spin resonancelike 
mode with E g — 2E S ~ 5 J with E s — (E m ) m - m denoting the lowest level in the spinon spectrum in terms of H s . 



7. Emergence of the nodal quasiparticle 



So far we have surveyed some novel properties of the SC state, which generally are non-BCS like. They involve 
either "high energy" (e.g., spin excitations above the spin gap E g ) or short-distance (e.g., the vortex core at a 
length scale ~ £) physics. However, in a sufficiently long wavelength and low energy regime, a typical d-wave 
BCS superconductor will eventually emerge, 29 where the physical properties will be essentially dominated by the 
conventional nodal quasiparticle excitations. 
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FIG. 11: A schematic nodal quasiparticle dispersion in the SC phase. The lineshape below the characteristic spinon energy 
E s = Eg/2 is very sharp as energetically the quasiparticle cannot decay into a spinon and a holon in the condensate. Above 
E a , however, such a decay is energetically allowed locally with the spinon and holon remain loosely confined at large distance. 
[From Ref. 31 ] 



A quasiparticle excitation is created by the electron c-operator. 29 According to the definition (17), it is composed 
of a pair of holon and spinon together with a phase vortex factor e ie as a "bound state" . Such a low-lying excitation 
is a "collective mode" , 29 which will be independent of the "mean-field" phase string model H stTing below the spinon 
gap E s = Eg/2. Besides the spinon gap in H Sl the gapless "phonon mode" for the condensed holons in Hh will be 
turned into the plasma mode by coupling to the external electromagnetic field — the Anderson-Higgs mechanism. So 
a low-lying quasiparticle should remain coherent since there are neither other low-lying states for it to decay into nor 
to scatter with. 

To further see how a quasiparticle is compatible with -ff s tring, one may examine the behavior of its three components. 
It has been seen that an isolated spinon excitation will cost a logarithmically divergent energy via the gauge field Afj . 
Similarly a single localized holon excited from the condensate will also create a divergent energy due to the global 
spinon response via the gauge field in H s . 23 Furthermore, a divergent energy is also associated with the creation 

of vortices by e* e . 23 However, as the bound state of these three fields, a quasiparticle excitation will no longer invite 
any nonlocal response from -ffstring neither from the holons nor the spinons. 23 ' 29 In other words, the holon and spinon 
are confined inside a quasiparticle, making it a stable "local" excitation with a finite excitation energy. The confining 
potential is logarithmic. 

The existence of the low-lying quasiparticle excitation may be also understood based on the ground state (87), 16 
which can be considered as the Bose condensation of the bosonic field D on the RVB background |RVB) with (^D^ ^ 0. 

Then a low-lying excitation may be constructed based on D\ x I'g)sCj with some smooth change in the D. Generally 
with achieving the phase coherence (84), one has 

£>|* G ) SC ~ ^crc^c^l^sc (109) 

ij a 

with 




(110) 



according to the discussion in Sec. 4.1.1. Therefore, such type of low- lying excitations, constructed with a smooth 
change in g^j, can always be described in terms of a pair of fermionic electronic (quasiparticle) excitations. In 
particular, due to the d-wave nature of fjij, the separation of the quasiparticlcs along the nodal line may be infinitely 
large. As "collective modes", the quasiparticles are not directly described by -Hiring, consistent with the fact that 
^string determines tp h and |RVB) in Eq. (75), but not the spinon backflow wave function Z in n. 

The effective Hamiltonian governing the motion of a quasiparticle excitation can be derived based on the original t-J 
model. 29,30 The renomalization effect from interacting with the background electrons will be represented by the mean- 
field order parameter A? c , etc. An equation of motion description of such a quasiparticle has been developed, 29 ' 30 
which shows the low-energy part of the spectral function is similar to a sharp d-wave BCS one, but the high-energy 
part (> E s ) is nontrivial as the composite (spin-charge separation) feature will show up, where the quasiparticle can 
decay into a pair of spinon and holon locally without costing much from the logarithmic confining potential. Figure 
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FIG. 12: Left: The single-electron spectral function in the one-hole case. The small bars mark the spinon spectrum -Ek+k "Mo- 
Right: The "quasiparticle" spectrum determined by the ARPES (open square) and by the spinon spectrum -Ek+k (solid curve). 
(From Ref. 9 ) 



11 illustrates the d-wave nodal quasiparticle dispersion schematically, where the quasiparticle spectral function has a 
sharp lineshape at 

-^qasiparticlc ^ -^spinon "I" -^holon (H-0 

Namely, the spinon gap E s (or spin gap E g ) will provide a protection for the coherent quasiparticles in the SC state. 
Indeed, in the half- filling limit, E s — > and the single-electron spectral function has only an incoherent "composite" 
part 9 as shown in Fig. 12. The more detailed results for the spectral function in the SC state will be presented in a 
forthcoming paper. 30 



B. Lower Pseudogap Phase (LPP) [Spontaneous Vortex Phase (SVP)] 

Due to the composite structure of the SC order parameter A sc [Eq. (94)], a regime may exist at T > T c where 
^ e «* ( r )^ = o but the Cooper pair amplitude A remains finite. Such a regime is known as the spontaneous vortex 
phase (SVP) or the lower pseudogap phase (LPP) of the phase string model. 26 ' 31 

The LPP is described by free spinon vortices, which are thermally excited and proliferate such that the phase $ s 
gets disordered. Namely, the LPP is an electron "fractionalized" state with the proliferation of unbinding spinons, 
and the main distinction between the SC phase and LPP lies in the phase (dc)coherence or spinon (de)confincmcnt. 

In the LPP, a finite A ensures that the spinon vortices are still well defined. In fact, ip h ^ means that the 
generalized GL equations [Sec. 4.1.2.] are also applicable in the LPP. Since the holon condensation persists in the 
LPP, the spin dynamics at E > E g should also remain qualitatively similar to the SC phase [Sec. 4.1.3.]. 



1. Phase diagram for the LPP 

A is composed of the RVB pairing A s and holon condensate ip* h , which disappear at some T and T v , respectively. 
We will see generally T v < T , as ip h ^ is always underpinned by the spin singlet pairing. Thus T v will represent 
the characteristic temperature for the LPP, whereas To defines the boundary of the so-called upper pseudogap phase 
(UPP) to be discussed in the next section. 

Note that without A s , the holon system would be a simple 2D hard-core boson problem according to Eq. (51) or (90), 

with the Kosterlitz-Thouless (KT) transition temperature for the holon condensation given by Tkt = 7r5 (2a 2 m/j) 
as shown in Fig. 13 by a dotted line (with th = 3J). 

However, the frustration effect of A s on the holon condensation will play a crucial role here. The spinon-vortex 
density is determined by n v — ^2 ma (llnolmo) /N ■ Due to the opening up of a spin gap E g in the holon condensation 
phase, n v is exponentially small for T <C E g . With increasing temperature, n v will monotonically increase until 
reaching the maximal number n™ ax — 1 — 5 at T = T where all the RVB pairs break up. 
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FIG. 13: Left panel: The characteristic temperature and magnetic field scales which determine the phase diagram of the upper 
and lower pseudogap and superconducting phases based on the phase string model. Here To and Hp G = Hpg(T = 0) for the 
UPP, T v and H® = H V (T = 0) for the LPP, and T c for the SC phase. Right panel: The phase diagram of the pseudogap 
regimes in the three-dimensional space of magnetic field, doping, and temperature. [From Ref. 31 ] 




FIG. 14: Two holons (open circles) exchanging positions can pick up a minus sign if an unpaired spinon (grey circle) is enclosed 
in between, which carries a ±7r flux tube. So the phase coherence of bosonic holons can be effectively destroyed to prevent the 
Bose condensation if there is a sufficient number of free spinon excitations in the spin background. 



Since each free spinon carries a ir fluxoid as perceived by the holons, the quantum phase coherence among bosonic 
holons will be violently interrupted if on average there is an excited (unpaired) spinon sitting between two neighboring 
holons (as illustrated by Fig. 14), where an exchange between a pair of holons can gain a minus sign in the wave 
function. In other words, the holon condensation must break down when the vortex density n v is comparable to holon 
density 6, far less than n™ ax at low doping. Such a consideration provides an estimate of the upper bound for T v as 31 

n v = 5. (112) 

Equation (112) can be also understood based on the "core touching" picture of spontaneously excited spinon vortices. 
Note that the average distance between excited spinons may be defined by l s = 2a/ ' ^J-Kn v . Since the characteristic 
core size of a spinon vortex is a c , then one expects that the "supercurrents" carried by spinon- vortices are totally 
destroyed when the sample is fully packed by the "cores" with l s = 2a c which also results in Eq. (112). 

The numerical result of the characteristic temperature T v in terms of Eq. (112) is plotted in the left panel of Fig. 
13 by the filled circles. It shows that T v is nested below Tkt at low doping and is dramatically reduced from Tkt 
with the further increasing of doping due to the frustrations induced by spinon excitations, which at larger doping 
remains nested below To of the UPP and eventually vanishes at srvb together with To (see the next section). 

An external magnetic field will break up more RVB pairs through the Zeeman effect (53) at a given temperature. 
By considering the Zeeman effect on the energy spectrum, E ma — E m — <j[i B B, the magnetic field dependence 
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FIG. 15: The magnetic field - temperature phase diagram of the pseudogap phases at doping S — 0.125. [From Ref. 31 ] 

of T v = T V (H) can be further obtained from Eq. (112) (note that we do not distinguish the magnetic field H and 
induction B as the magnetization is usually very weak here). Or conversely, for each T < T v (0) there is a characteristic 
field H V (T) at which the LPP phase is destroyed [Fig. 15]. H° = H V (T = 0) determined this way is shown in the left 
panel of Fig. 13 by the filled squares. 

For comparison, the SC temperature T c [Eq. (108)] is shown as the dashed curve in the left panel of Fig. 13. 
Furthermore, in the mixed state below T c , by including the Zeeman energy, E g is reduced to E* — E g (B) ~ 2fi B B such 
that one can estimate T C (B) by using a simple relation T C (B) ~ E*/4k&. Then in turn one may define H* 2 = B(T C ), 
which is shown in Fig. 15 by a dashed curve. Note that H* 2 so defined will vanish at T Cl resembling the conventional 
H C 2 in a BCS superconductor. However, since free spinon vortices are generally present at H > H* 2 , H* 2 is a crossover 
field which no longer has the same meaning as H C 2 in a conventional BCS superconductor. Roughly speaking, the 
Abrikosov magnetic vortices are expected to be present below H* 2 where the spontaneous spinon-vortices, generated 
by the Zeeman effect, are still loosely paired, whereas the vortex unbinding occurs above H* 2 . So H* 2 defines a 
crossover between two types of vortex regime. The numerical result shows that [i*b,H* 2 {0) ~ E g (B = 0, T = 0)/2, 
which results in H* 2 (0) ~ 3T C (Tesla/Kelvin) , according to Eq. (108). 



2. Nernst effect 



Base on the London- like equation (98), using the steady current condition 

a t j = o 

and the electric field E = —d t A e in the transverse gauge, one finds 26 

E = z x d> n (n+v-L - nlv-) 



(113) 



(114) 



where denotes the density of spinon vortices and antivortices with drifting velocity v± along a direction perpen- 
dicular to the electric field. As illustrated by Fig. 16(a), the electric field and the drifting of vortices and antivortices 
must be balanced according to Eq. (114) in order to avoid the system being accelerated indefinitely with d t J ^ 0. So 
the applied electric field will drive the vortices and antivortices moving along a perpendicular direction with opposite 
velocities: v + = — v if the vortices and antivortices are not polarized by the external magnetic field, i.e., n+ = 

The Nernst signal ejv is the electric field measured along the y-direction when spinon vortices and antivortices are 
both driven by a temperature gradient in the same direction along the x-direction: 



E,, 



e N 



-V X T 



(115) 



Such a case is shown in Fig. 16(b), where the spinon-vortices and -antivortices move along the i-direction with 
v + = v_ = v. To have a finite E in terms of Eq. (114), then the vortex density has to be polarized by the 
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FIG. 16: Schematic picture of (a) the flux flowing under an electric field in the LPP, which can lead to a flux flow resistivity as 
well as the spin Hall effect, and (b) the flux flowing under a temperature gradient VT, which must be balanced by an electric 
field E and thus leads to the Nernst effect. 



external magnetic field B —Bi according to the "neutrality" condition B = <fi Q (rt+ — n v ) such that 

E = Bxv. (116) 

Suppose S0 is the transport entropy carried by a spinon vortex and rj s is its viscosity such that the drift velocity 
v s can be decided by s^VT = —iJs v - Then one has e^r = B—. On the other hand, in the absence of the temperature 
gradient, a charge current can also drive a transverse motion of spinon vortices and antivortices along opposite 
directions, i.e., v± = ±v, such that an electric field is generated along the current direction according to Eq. (114), 
leading to a finite resistivity due to the presence of free vortices, which is given by 

P=-<t>l (117) 

Vs 

with n v = + n~ . This formula is familiar in the vortex flow regime of a conventional superconductor. Then, by 
eliminating rj s , one obtains 31 

a xy = — = -j^. (118) 

P <Po n v 



What really makes the Nernst transport unique in the present theory is that the transport entropy here is 
associated with the spin degree of freedom due to its free 5 = 1/2 moment, instead of a normal core in a conventional 
BCS superconductor, which is given by s$ = /cb {In [2 cosh ((3/i B B)] — P^B tanh (f3fi B B)} . The temperature and 
magnetic field dependence of a xy is shown in Fig. 17. 31 The magnitude of such a quantity is comparable to the 
experimental data, implying that the transport entropy due to the free moment in a spinon-vortex is capable to 
produce the Nernst signal observed experimentally. 32 



3. Spin Hall effect 



A unique prediction related to the spinon-vortex motion driven by an external electric field is the existence of a 
conserved dissipationless spin Hall current. 19 As shown in Fig. 16(a), vortices can be driven by an in-plane electric 
field to move along the transverse direction. Since each vortex carries a free moment, if these moments are partially 
polarized, say, by an external magnetic field via the Zeeman effect, then a spin Hall current can be generated along 
the vortex motion direction. The spin Hall conductivity is determined as follows: 33 
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FIG. 17: The quantity a xy (j> /d = eN<f> / pd, which is related to the Nernst signal without involving the viscosity coefficient, is 
shown as a function of the magnetic field and temperature (the inset) at a given doping. Here d — 7.7 A is the distance between 
two Cu02 layers. [From Ref. 31 ] 



r* = ( B 

" 9Vb \ih,<j> 



(119) 



which only depends on the intrinsic properties of the system like the uniform spin susceptibility x s , with the electron 
^-factor g ~ 2. It is important to note that the external magnetic field B applied perpendicular to the 2D plane reduces 
the spin rotational symmetry of the system to the conservation of the S z component only, satisfying + V • 3 s =0. 
Thus the polarized spin current J s is still conserved and remains dissipationless as the current of its carriers - vortices 
is dissipationless in the LPP. In contrast, the charge current remains dissipative as p ^ 0. 



4- Magnetization 

Inspired by the experiment, 32 the diamagnetism has been also studied. The total magnetization can be expressed 
as 34 

Mtot = A/dia + M para (120) 

in which Mdi a and Af para stand for the orbital diamagnetism from the vortices and the paramagnetism from the 
Zeeman coupling, respectively. Based on the mutual Chern-Simons formulation of the phase string model outlined in 
Sec. 3.2, with a systematic description of multi-spinon excitations and interaction between them, the nonlinear effect 
of the magnetization vs. magnetic field can be effectively treated. 34 

The magnetic field and temperature dependence of the total magnetization at different doping concentrations as 
well as the diamagnetism part Mdj a at S = 0.125 are shown in Fig. 18 based on a mean-field approximation in the 
mutual Chern-Simons description. 34 Note that M para = \ S B with x s as the uniform spin susceptibility to be discussed 
later in the weak field limit. 



5. Magneto-resistance 



Due to the deconfinement of spinons, the quasiparticles are no longer the stable low-lying excitations in the LPP. 
The origin of the dissipation comes mainly from the flow of the spinon vortices. In fact, the resistivity (117) is similar 
to the flux-flow resistivity in a Type II superconductor except that n v in general is not simply proportional to the 
external magnetic field B. Namely, the spinon vortices can be spontaneously (thermally) generated with n v ^ 0, such 
that p ^ even at B = 0. The resistivity p{B) can be then expanded as 34 



p(B)= P (0) [l+ 7 S 2 + 0(i? 4 )] 



(121) 
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FIG. 18: The B dependence of the total magnetization M = Mdi a + Af para at various doping concentrations: (a) S — 0.078, (b) 
5 = 0.125, (c) S = 0.188; (d) the diamagnetism M dia at 8 = 0.125. [From Ref. 34 ] 
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FIG. 19: The magneto-resistance coefficient 7 vs temperature for the magnetic field B which is either perpendicular or transverse 
to the ab plane. [From Ref. 34 ] 



where the odd power terms of B vanish due to the symmetry p(B) = p(—B). Suppose that the dependence of the 
viscosity rj s on B is negligible, then the quadratic coefficient 7 can be expressed as 

_ p(B)-p(0) _ n„(£?)-n„(0) 
7 p(0)B 2 ~ n v (0)B 2 ' 1 ' 

The coefficient 7^ and 711 , with the external magnetic field B perpendicular and parallel to the 2D plane, respectively, 
can be calculated numerically as shown in Fig. 19. An important prediction of the present theory, as shown by 
Fig. 19, is that 711 is comparable to 7j_. This is a rather unusual case for a vortex-flow-induced resistivity, since 
normally the in-plane vortices are always created by the perpendicular magnetic field in a Type II superconductor, 
where the vortex-flow-induced resistivity only exhibits field-dependent magneto-resistivity for the component of B 
which is perpendicular to the plane. But in the present theory, vortices are tied to the free spinons. Since the 
latter can be created by the Zeeman term with the external magnetic field pointing at any directions, the former 
can thus be created by the in-plane field as well. The present mean-field-type treatment of p may not be expected 
to be quantitatively accurate in view of possible corrections from the fluctuations, but the existence of an anomalous 
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FIG. 20: The characteristic temperature To of the UPP versus 8/xnvB- Solid line: the present theory; Full squares: determined 
from the uniform spin susceptibility x s m LSCO compound; 36 Open circles: determined from the c-axis magneto-resistivity 
(p c ) measurement in Bi-2212 compound; 37 The dashed line shows the scaling relation of the zero-temperature critical field Hp G 
with To as predicted by the theory. Inset: the critical field Hpg as a function of temperature at 5 = 0.22. The experiment data 
from the c-axis transport in Bi-2212 (Ref. 37 ) are also shown by the open and full squares. [From Ref. 22 ] 



transverse magneto-resistivity with — y || comparable to 7j_ remains a very peculiar prediction: 



AtT>T v , the bosonic holons will be in a non-degenerate regime where the quantum coherence is destroyed by the 
excited spinons via A s . By contrast, the bosonic spinons can still maintain their quantum coherence up to To where 
the spin singlet pairing order parameter A s eventually vanishes. The UPP as characterized by A s ^ is defined at 
T v < T < T , whose key features will be described in the following. 



Figure 20 shows To (solid curve) as a function of doping with J = 1350 K. The experimental data determined by the 
uniform spin susceptibility measurement in LSCO 35,36 (see the discussion in the next section) are shown by the full 
squares. Furthermore, the open circles are independently determined from the c-axis transport 37 in the overdoped 
regime. Note that here J g = (1 — 2gS) vanishes at S = xrvb = 1/ (2<?) and the curve T versus S/xyo/b in Fig. 20 is 
not sensitive to the choice of g. Based on the above experimental data 35-37 xrvb is fixed at 0.25 such that g = 2. 

Due to the bosonic RVB origin of the UPP, the Zeeman effect of an external magnetic field can effectively destroy 
the singlet pairing of spins in the strong field limit, which is the only direct field effect on the RVB background. 

With incorporating the Zeeman term (53), one can obtain the zero-temperature "critical" field Hp G at which A s 
vanishes: 



C. Upper Pseudogap Phase (UPP) 



1. Phase diagram for the UPP 



Based on the self-consistent solution of H s , the characteristic temperature Tq at which A s — > is given by ; 
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FIG. 21: (a) The calculated uniform spin susceptibility \s scaled with the maximum x™ ax a ^ To versus T/To, which follows 
an approximately doping-dependent curve, (b) The theoretical \ s &t half-filling (solid) and the one obtained by the high 
temperature series expansion (HTSE). The latter fits the experimental scaling curves 35 ' 36 very well. [From Ref. 22 ] 



In Fig. 20, ln3/j, B Hp G / hs is plotted as the dashed curve which scales with the zero-field To fairly well, which means 
H B Hp G ~ l.lfceTo. The temperature dependence of the "critical" field Hpq(T) is shown in the inset of Fig. 20 at 
S = 0.22, together with the experimental data obtained from the c-axis magneto transport measurements. 37 We see 
that the high-temperature experimental data (open circles) fit the theoretical curve very well without any additional 
adjustable parameter. Furthermore the zero-field To determined by the same experiments is also in good agreement 
with the theory as shown (open circles) in the main panel of Fig. 20. But one may notice that the experimental 
Hpq(T) starts to deviate from the theoretical curve in the inset (full squares) as the temperature is further lowered 
and saturated to approximately the half of the predicted number (which implies Mb-^pg — /cbTo/2). However, such 
a deviation occurs only for those data (full squares) which have been obtained by extrapolation in the experimental 
measurement 37 and therefore may not be as reliable as the higher temperature ones (open squares) in the inset of 
Fig. 20. 

Finally, the three-dimensional phase diagram of the UPP, together with the LPP and SC phase, in the parameter 
of magnetic field, doping concentration, and the temperature is summarized in the right panel of Fig. 13. Hpq vs. T 
at 6 = 0.125 is also plotted in Fig. 15. 



2. Uniform spin susceptibility 

The spin singlet pairing (RVB) nature of the UPP is clearly manifested in the uniform spin susceptibility \ s given 
in the main panel of Fig. 21(a) at different doping concentrations. 27 Note that x s reaches a maximum value xT ax 
at temperature T where the RVB order parameter A s vanishes. At T > T , x s follows a Curie-l/T behavior as 
spins become free moments at the mean-field level. The curves in Fig. 21(a) are presented as X s /x™ ax versus T/Tq, 
which approximately collapse onto a single curve independent of doping. The comparison with experiment has been 
discussed in Ref. 22 . 

In Fig. 21(b), the calculated x s versus T at S = is shown together with the high temperature series expansion 
(HTSE) result 38 . It is noted that the experimental scaling curve actually coincides with the half- filling HTSE very 
well. 35,36 Thus one can clearly see the overall qualitative agreement between the bosonic RVB theory and the experi- 
ment from Figs. 21(a) and (b). Note that the mean-field x s deviates from the HTSE result prominently around T 
where the latter is a much smoother function of T. It reflects the fact that To is only a crossover temperature and the 
vanishing A s does not represent a true phase transition. Obviously, the amplitude fluctuations beyond the mean-field 
A s have to be considered in order to better describe x s in this regime. T determined in the mean-field theory is quite 
close to the HTSE result, indicating the crossover temperature itself can still be reasonably decided by the mean- field 
bosonic RVB description given above. 
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FIG. 22: Left panel: 17 Ti/ 63 Ti vs. temperature at different doping concentrations in the upper pseudogap phase of the b-RVB 
state. The dashed line shows the result of an f-RVB state (tt flux phase) at half-filling. The inset shows the non-Korringa 
behavior of l/ 63 TiT in the b-RVB state at various dopings. Right panel: I/T2G vs. temperature in the upper pseudogap phase 
below To at different doping concentrations. [From Ref. 22 ] 



3. Spin-lattice relaxation and spin-echo decay rates 

The NMR spin-lattice relaxation rate of nuclear spins is determined by 22 

1 _ 2k B T v 

where the form factor T(q) 2 comes from the hypcrfinc coupling between nuclear spin and spin fluctuations. Due 
to the fact that the T(q) 2 for planar 17 0(2) nuclear spins vanishes at the AF wave vector Qaf> while for planar 
63 Cu(2) nuclear spins is peaked at Qaf, a combined measurement of l/ 63 Ti and l/ 17 Ti will thus provide important 
information about the AF correlations at low frequency u — > 0. 

Based on the mean-field equation, 22 the calculated spin-lattice relaxation rates, l/ 63 Ti and l/ 17 Ti, for the planar 
copper and oxygen nuclear spins are presented in the left panel of Fig. 22. It shows that the ratio T\j Ti, which 
is a constant above To, starts to increase with reducing temperature below To. At lower temperature, T/To < 0.5, 
such a ratio arises sharply. For example, 17 Ti/ 63 Ti diverges at 6 — as a true AFLRO exists at T = 0; And it can 
still reach about 100 in the low temperature limit at S = 0.125, all qualitatively consistent with the experimental 
observations. As pointed out above, such behavior clearly demonstrates that strong low-lying AF correlations around 
Qaf develop in the UPP, leading to the simultaneous enhancement of 1/ 63 Tl and the cancellation in l/ 17 Ti. In the 
inset of the left panel in Fig. 22, l/ 63 TiT has been plotted, which is also qualitatively consistent with the experiment, 
but deviates from the conventional Korringa behavior l/ 63 TiT ~ const for a Fermi liquid system. By contrast, the 
ratio 17 Ti/ 63 Ti in an f-RVB mean-field state (the n flux phase) at half-filling remains flat over the whole temperature 
region as shown by the dashed line in Fig. 22, indicating the absence of any significant AF correlations around Qaf 
in the pseudogap regime of the fcrmionic RVB state. 

The spin-echo decay rate I/Tzg, which is related to the static AF correlations via the real part of spin susceptibility 
function, is also examined in the right panel of Fig. 22. It shows that I/Tzg begins to increase with reducing 
temperature below T>. Such behavior has been also observed in the experiment, which once again clearly supports 
the picture that the strong AF correlations start to develop in the UPP. 

By comparison, the corresponding magnetic properties in the LPP are presented Fig. 23. In the main panel, the 
uniform spin susceptibility shows a true "spin gap" behavior, in contrast to the "scaling" curve shown in the UPP in 
Fig. 21 where x s m the doped regime roughly behaves like that at half-filling — in the latter case x s saturates to a 
constant at T = 0. In the LPP, these x s ' s can drop below that at d = and vanish at T = 0. Furthermore, l/ 63 Ti 
also decreases with temperature (see the left inset of Fig. 23), as opposed to the behavior in the UPP, indicating the 
appearance of the spin gap over whole momenta. On the other hand, although the low-energy spin fluctuations are 
gapped, the static AF spin-spin correlations as described by the real part of spin susceptibility function still remain, 
as reflected by I/T2G shown in the right inset of Fig. 23, where the monotonic increase of I/T2G in the UPP (Fig. 
22) is replaced by the saturation at the LPP. 
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FIG. 23: Uniform spin susceptibility in the lower pseudogap phase at different dopings including half-filling. The left insert 
shows l/ 63 Ti and the right insert I/T2G with the same symbols as in the main panel. [From Ref. 22 ] 



D. Incoherent "Normal State": Classical Regime 



As shown above, at T > T , the mean- field RVB order parameter A s vanishes and to the leading order of ap- 
proximation the spins are localized at the lattice sites with very weak correlations between them. The residual AF 

2 

superexchange coupling will come from Af - which has been neglected in the minimal phase string model. 

In this regime, the charge dynamics will be highly nontrivial as governed by Hh- According to Eq. (29), an isolated 
spinon excitation will behave like a 7r-fiux tube as perceived by the holons and thus provides a strong, unconventional 
charge scattering source. So at high temperature where a lot of spinons (<~ 1 — 5) are thermally excited, one expects 
a severe intrinsic frustration effect exerted from A s . on the holons. 



1. Novel scattering mechanism for the charge carriers 

To see how the spin dynamics influences the charge degree of freedom via A s , one may write down the propagator 
for A s [Eq. (91)] 39 

r0 

D^(d,iu n ) = / dre— (T T A«(q,r)^(-q,0)) 
Jo 

= -(^-^)^X"(q,i«n) (126) 

Define the local flux per plaquette (surrounding a lattice site) $□ = a 2 i- (V x A s ) . Its total strength is generally 
determined by 

((<I>&) 2 ) = f ^l]T47r 2 ^(q,c) (127) 
J q 

where the spin structure factor S zz (q, u>) = tt^ 1 [1 + Imx zz (q,w). 

In particular, at T > To, there is no more significant AF correlations among spins as A s = and one finds 

S zz (q,uj) = ^^/ {3 ~ 6) ^~^ S{uj). Thus the corresponding gauge flux fluctuation becomes truly static with the weight 

~ -k\J ' ( 3 ~^^ 1 ~^) concentrating at w = 0. At low doping, (($^) 2 ) is comparable to the simple picture that each 

excited spinon contribute to a flux of order of n, which represents the maximal frustration effect that the holons can 
experience in the phase string model. 39 
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FIG. 24: Solid curve: the holon density of states (DOS) determined by Hh with the flux strength |$ n = (1 — 5)ir at 5 = 0.125 
and the dash-dotted curve: with the reduced flux strength N& n = 0.47T. By comparison, the dashed curve represents the 
flux-free limit, while the dotted curve corresponds to the case in the presence of uniform ir flux per plaquette. [From Ref. 39 ] 




FIG. 25: Left: The real part of the optical conductivity at 5 = 0.0625 with |$p| = (1 — 0.625)7r. Right: The optical conductivity 
at S = 0.125 with |<E> n | chosen from |$ n | max = 0.875tt to 0.2tt at a fixed T = 0.5t h . [From Ref. 39 ] 



Figure 24 illustrates 39 how the holon density of states (DOS) will get reshaped by A s or $p. The quenched method 
is used to average over the static random flux configurations for $g = ± |<f>n| on the lattice plaquette. It shows 
how the DOS drastically reshaped by the gauge field, i.e., the suppression in the high-energy (mid-band) DOS, as 
compared to the flux free case (dashed curve). Note that the dotted curve in Fig. 24 represents the DOS for the case 
of a uniform tt flux per plaquette, which looks similar to the random flux case |d>rj| = 0.875 except that the momenta 
remains well defined in a reduced Brillouin zone in contrast to a strong mixing of momenta over a wide range by the 
scattering effect in the latter. 



2. Optical conductivity 

Figure 25 shows the real part of the calculated optical conductivity. 39 The main feature of the spectral curves at 
various temperatures T > To is that there is generally a two-component structure with a usual low-energy Drude 
component (~ l/^ 2 ) and a mid-infrared resonancelike peak around the energy scale w^ir ~ 4t/j. Furthermore, such a 
mid-infrared peak will actually smoothly evolve into the 1/uj behavior with reducing |$ n |, as clearly illustrated in the 
right panel of Fig. 25 at a fixed holon concentration <5 = 0.125 where the mid-infrared resonancelike peak at smaller 
<& p |'s becomes softened and finally behaves like al/w tail in the regime ~ 2th — &th with the weight shifting towards 
the lower energy. 

The origin of the mid-infrared resonance has been one of the most intriguing optical properties in the underdoped 
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FIG. 26: Optical conductivity at different temperatures below To. A new lower mid-infrared peak emerging at u^ir which is 
weighted between the magnetic energy scales E g and _E" ppcr . [From Ref. 39 ] 



cuprates. Normally a photon with the momentum q ~ cannot excite a high-energy particle-hole pair involving a 
large momentum transfer due to the momentum conservation law, unless there is a scattering mechanism to strongly 
and widely smear the momentum. This is difficult to realize in a conventional elcctron-collcctive-modc coupling 
mechanism. The phase string model provides an alternative scattering mechanism due to the strong correlation effect. 

We have already seen that the effect of A s results in a double-peak structure in the holon DOS (Fig. 24). In 
contrast to the uniform n flux case shown in the same figure, which also has a double-peak structure, the high-energy 
inter-peak transition at q — > becomes possible in the random flux case due to the mixing between the small and 
large momenta by the strong scattering via A s . This is the origin for the mid-infrared peak found in Fig. 25. 

Finally, as a comparison, the low-energy optical conductivity at T < T can be obtained by using the perturbative 
method. 39 Such an approach is meaningful in the regime where the spin fluctuations are substantially suppressed, 
which in turn results in the weak fluctuations of A s according to Eq. (126). As shown in Fig. 26, a prominent 
suppression of a'(uj) at low-w is present at T = with a second "mid-infrared" peak emerging around w^ ir <~ 0.75 J 
which sits somewhat between the two characteristic magnetic energy scales, E g and E^ ppcl , as marked in the figure. 
Note that such a new energy scale in the low-w optical conductivity merely reflects some weighted energy scale based 
on the magnetic lm\ zz - With the increase of temperature, the "gap" at low energy in u'(lo) is quickly filled up by 
the thermal excitations as shown in Fig. 26. The lower "mid-infrared" peak feature remains around cj^ ir at low 
temperature throughout the LPP below T v . Note that T v is between T c and T , and the dashed curve at T = 0.4T is 
obtained by supposing that T > T v where lmx zz behaves differently. 22 As compared to the solid curve at the same 
T = 0.4Xo, which corresponds to the case inside the LPP, the overall difference is small except for the vanishing the 
lower "mid- infrared" peak [Fig. 26]. 



3. Density- density correlation function 

The mid-infrared resonance peak of the q = optical conductivity has been attributed to a large-w transition 
between the double peaks of the holon DOS [Fig. 24]. In the following we discuss an independent probe of such a 
peculiar DOS structure by studying the density-density function at finite momentum q and energy ui, and compare 
the results with the exact numerical calculations. 

The imaginary part of the (retarded) density-density correlation function Cd(q,u>) is presented in Fig. 27 (solid 
curves), which evolves distinctively with different momenta. Note that Cd(q, w)/5tt is shown in the figure because it 
is quantity roughly doping independent. For comparison, the exact diagonalization results 40 are presented as dotted 
curves. It is interesting to see that the overall w-peak feature of the calculated density-density correlation function 
is in qualitative and systematic agreement with the numerical one at different q's without fitting parameters (here t 
is simply set at th as the mid-infrared feature peaks around <~ At in the numerical calculation). Such a consistency 
between the present effective theory and the exact diagonalization provides another strong evidence that the gauge- 
coupling boson model (51) correctly captures the high-energy charge excitations in the t-J model and large-C/ Hubbard 
model. 
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FIG. 27: The calculated density-density correlation function Cd(q,cj)/ 'Stt at different momenta (solid curves) with T = 0.5th- 
The exact diagonalization results 40 are shown as dashed curves. [From Ref. 39 ] 



4- Scattering rate 
Experimentally the scattering rate is normally defined by 

-h - (t- ) Re 

t{uj) \ 4ir I 

which is determined by the measured optical conductivity. Here uo p denotes the plasma frequency, which in the present 
case is given by ui p — ^J%Tte 2 8th. 

In the left panel of Fig. 28, 1/t(u) based on the calculated <r(w) is plotted as a function of u) in different temperatures 
at 6 = 0.125, with |$ n | = 0.47T which corresponds to the case where the high-u> optical conductivity looks more like a 
1/lu behavior (Fig. 25). Here one finds that 1/t(u) increases monotonically with lo and is roughly lincar-u> dependent 
over a wide us region at w > fee^o- Note that generally the w-dependence of 1/t(uj) at higher energies is closely 
correlated with the evolution of the aforementioned mid-infrared feature, as shown in the inset of the left panel of 
Fig. 28. 

In particular, one sees a parallel shift of 1/t{lo) with increasing temperature at low-w, which implies a linear- 
temperature dependence of the dc scattering rate. The dc scattering rate I/t^c can be determined by extrapolating 
\/t{lo) to lo = 0. The obtained dc resistivity based on the Drude formula p dc = (uip/in) t7 c = l/er'(0) is shown in 
the right panel of Fig. 28 which is indeed quite linear over a very wide range of temperature at T > Tq. 

It is important to note that cr'(O) oc j3 over a very wide range of the temperature at T > T where the Bosc 
distribution factor n(£ m ) <C 1, i.e., in the classical regime of the bosons (£ m is the holon energy spectrum). The 
corresponding scattering rate H/t^c ~ 0.7fcgT for the case shown in the main right panel of Fig. 28, whose slope is 
slightly |$ n | dependent as indicated in the inset. Indeed, as discussed above, the bosonic degenerate regime for the 



a{u) 



(128) 
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FIG. 28: Left: The scattering rate 1/t(cj) denned by Eq. (128) at various temperatures between To ~ 0.25th and 2Tb which 
show a rough linear-aj dependence over a wide range at uj > feTo. Inset: 1/t(lu) vs. u> at different |$rj|'s corresponding to Fig. 
25. Right: The dc resistivity p dc = l/a'{u)\ U r^o as a function of temperature which is fit by a straight line showing the linear-T 
dependence. Inset: p dc at different |$rj|'s which all show good linear-T behavior with slightly different slopes. [From Ref. 39 ] 



holons already ends up at T v , i.e., at the boundary of the LPP/SVP (Fig. 13). At T > Tq, totally 1 — 5 randomly 
distributed 7r-flux tubes are perceived by the S holons and the latter behave like classical particles. One expects this 
anomalous transport be smoothly connected to the Brinkman-Rice retracing path regime 41 in the large T limit. 

The dc scattering rate fi/ T dc ~ 2fc^T has been previously obtained 42 by the quantum Monte Carlo numerical 
method, where the starting model is a system of interacting bosons coupled with strong Gaussian fluctuations of the 

static gauge field of the strength ^(^q) 2 ^) • Note that (^(^n) 2 ) used in the Monte Carlo simulation 42 is about the 

same order of magnitude as in the above case and, in particular, it is temperature independent in contrast to a linear-T 

dependent ^(^n) 2 ^) predicted in the slave-boson U(l) gauge theory 18 which was the original motivation for such a 

Monte Carlo study. 



E. Low-Doping AF State: Beyond the Minimal Model 



At half-filling, the antiferromagnetism can be well described by H s . However, even in the presence of a very dilute 
hole concentration, an AFLRO would be immediately destroyed due to the opening up of a spin gap E g cx 5 J as 
predicted by the minimal phase string model, as illustrated by the phase diagram in Fig. 13. In the following we 
discuss a modified phase diagram by taking into account of a new topological excitation in this regime. 

The motion of holes will generally induce the irreparable phase string effect (Sec. 2.2.). In the dilute limit of the 
hole concentration, the phase string effect should mainly influence the hole dynamics, without drastically affecting 
the spin part which is AFLRO ordered in the ground state. It turns out that the holes can be self-localized by the 
phase string effect here. 9 ' 43,44 Without the condensation of the holons, then a spin gap E g cx S J will no longer exist 
in this dilute doping regime. 

Mathematically, a Z2 topological excitation 45 (meron) is allowed by the phase string model 43,46 

b ia -> b la ef^ (129) 



h\ - hie*™ (130) 

where -d^ — ±Imln(zi — zQ) with the core position z° either inside a plaquette or on a lattice site. Such a meron 
can be "nucleated" from the AF state where the spinons are Bose condensed with (bi a ) 7^ and the holon becomes 
a topological vortex with a logarithmically divergent energy. 43,46 A holon must be then forced to be "confined" to a 
meron to form a truly stable hole object, known as a hole dipole. 43 ' 46 Two typical hole-dipoles of minimal size are 
sketched in Fig. 29: It can be shown that the effective mass of the induced meron is infinity such that the hole-dipole 
object is self-trapped in space. 44,46 

In such a self-localization regime, the kinetic energy of the holes is suppressed. Without the balance from the 
kinetic energy, the low-energy physics in this regime will be determined by potential energies. The latter will then 
decide various competing orders in this low-doping insulating phase. The AFLRO should persist if a weak interlayer 
coupling is considered. 43 But the freedom in the directions of the hole-dipole moment will lead to the reduction of 
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FIG. 29: (a) Type A dipole with a holon and a spin sitting at two poles at nearest-neighboring lattice sites; (b) Type B dipole 
with one pole at the center of a plaquette. [From Ref. 46 ] 




FIG. 30: Phase diagram at low doping x: a dual confinement-deconfinement occurs at a quantum critical point x c ~ 0.043. 
The Neel temperature Tjv vanishes at xq ~ 0.03. Tf and T g denote characteristic spin freezing temperatures, and T c is the 
superconducting transition temperature. [From Ref. 43 ] 



the Neel temperature T/v as shown in Fig. 30. Based on the hole-dipole picture and the renormalization group (RG) 
calculation, 43 the critical doping xq ~ 0.03. Beyond xo or Tn, the system is in a cluster spin glass phase with the 
dipole moments being quenched randomly in space. With the further increase of doping, the sizes of hole dipoles will 
get larger and larger, and eventually a deconfincmcnt can occur at a critical doping 6 = x c ~ 0.043, beyond which 
single holons will be unbound from their anti-vortex partners and experience a Bose condensation leading to a finite 
E g cx VS-x c J. 43 

Some interesting properties including the thermopower, variable-range-hopping resistivity, and dielectric constant, 
etc., in this self- localization regime have been also discussed. 44 Furthermore, a possible stripe instability has been 
explored. Note that in the above it is assumed that the hole dipoles are self-trapped uniformly in space. However, if 
there is no impurities or disorders, the uniform distribution of the self-trapped hole-dipoles may not be stable against 
the formation of the stripes due to the long-range dipole-dipole interaction. 46 ' 47 

So in the low-doping regime, the minimal phase string model in Eqs. (51) and (52) can be generalized to 

H h = /,^T(<' 1 " )//>, • U. r. (131) 

m 

H s = -J. J2 blb)_ a + H.c. (132) 
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to include a Z 2 gauge field fi^-: 

E n « = { ± n 7r ' (133) 

which is allowed by the general construction of the phase string model and compatible with the bosonic RVB pairing 
A s . Normally the core energy of a Z 2 vortex is too big in the superconducting phase, but a Z 2 vortex excitation can 
become important at <5 < x c in the AF state. Furthermore, a quasipaticle excitation discussed in Sec. 4.1.7. may 
be equivalently considered as a bound state of a spinon, a holon, and a Z 2 vortex. The details will be presented 
elsewhere. 



V. SYNTHESIS AND PERSPECTIVES 

In this brief review, I have surveyed a systematic effort in the study of a doped antiferromagnet, which may be 
relevant to the mechanism of the high-T c cuprates. The core of this approach lies in the so-called phase string effect, 
which has been mathematically identified based on the t-J model. It is by nature a frustration effect induced by the 
motion of doped holes in an antifcrromagnetic spin background. 

Such a frustration effect on the spin degrees of freedom differs fundamentally from an ordinary frustrated antiferro- 
magnet in the presence of geometrically "mismatched" spin interactions, e.g., the next nearest neighbor superexchange 
coupling. The key distinction is that the frustration in the former is dynamically driven and mutual between the spin 
and charge degrees of freedom. Namely the extent that the spins get frustrated crucially depends on the charge behav- 
ior and vise versa. In different doping, temperature, magnetic field, or other parameter regimes, the spin and charge 
parts will then behave differently in order to minimize the total free energy. For example, in the dilute hole limit, the 
superexchange energy of spins will dominate and with maintaining longer range antifcrromagnetic correlations the 
kinetic energy of doped holes can get severely frustrated by the phase string effect, resulting in their self-localization 
at low temperature; At higher doping, to gain the kinetic energy of the doped holes, however, the spin correlations 
can be "forced" to become short-ranged via the phase string effect, and the spin background becomes a spin liquid 
state. The superconducting phase coherence and nodal quasiparticle excitation are protected by the spin gap of such 
a spin liquid state. 

The mathematical description of the phase string effect is rather simple, which is basically represented by a sequence 
of signs [Eq. (11)] picked up by the nearest neighboring hoppings of the holes in a Heisenbcrg spin background. It 
depends on the path of the hole hopping as well as the polarizations of those spins exchanged with the hole during its 
motion. It is thus both geometric and dynamic, which weights each motion path of the holes. We have seen that such 
a phase string is irreparable in the sense that the system cannot generate other signs to compensate it at each path. 
In fact, the Heisenberg superexchange interaction respects the Marshall signs, so a phase string as the disordered 
Marshall signs caused by hopping cannot be "self-healed" through the superexchange process. 

Such an irreparable phase string effect identified as the most important frustration effect in the doped antiferromag- 
net is singularly sensitive to any perturbative treatment, since Eq. (11) will change sign for a fluctuation with merely 
one additional or less J, spin exchanged with the hole on a given path, no matter how long the path is. Fortunately 
a unitary transformation exists in the t-J model which can precisely keep track of such a phase string because it 
essentially involves the counting of the exchanges occurring between the holes and spins during the their travelling. 
Then in the new representation after the unitary transformation, known as the phase string formalism, the t-J model 
presumably becomes less singular and perturbatively treatable. 

The t-J model in the exact phase string formalism is a topological gauge model, in which the phase string effect 
is precisely described by a pair of mutual Chern-Simons gauge fields in two dimensions. Without these gauge fields, 
the model reduces to a full bosonic one free from any "sign problem" . In other words, the nontrivial spin and charge 
dynamics will be governed by the topological gauge fields which precisely reflect the phase string effect. Thus, the 
exact phase string formalism of the t-J model provides a unique starting point to study the doped antiferromagnet. 

Based on the precise topological gauge structure and the good understanding of the half-filling phase, an effective 
minimal phase string model working for small doping can be then constructed as given by Eqs. (51) and (52). Despite 
its novel looking, this is a rather simple model where two bosonic matter fields, spinless holons and neutral spinons, 
interact with each other by perceiving the opposite species as a 7r flux fluxoid, and the spinons form the RVB pairing 
whose amplitude is self-consistently determined. 

Figure 31 summaries the global phase diagram of this minimal model. The "normal state" at T > To is an 
incoherent "classical" regime for both spinons and holons, where the spinons are weakly correlated, while the diffusive 
holons are maximally scattered by the ±7r flux fluxoids bound to the randomly distributed spinons, leading to a 
linear-temperature scattering rate. At T < T , the spinons start to form the singlet RVB pairing, and the short-range 
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FIG. 31: The global phase diagram of the minimal phase string model. 



antifcrromagnetic correlations become enhanced with the reducing temperature, as clearly indicated by the NMR spin 
relaxation rate and spin-echo decay rate. This regime is known as the upper pseudogap phase, which is superexchange 
energy driven and continuously evolves into the AFLRO state at half- filling and T = 0. The holons remain diffusive in 
this phase, although the scattering from the gauge field gets reduced with decreasing temperature as more and more 
spinons are RVB paired with the cancellation of their ±7r fluxoids. Eventually at T < T v , the gauge flux is reduced 
so much that the bosonic coherence between the holons can be finally established, which defines the lower pseudogap 
phase that is obviously kinetic energy driven. In this phase, with the holon condensation, a spin gap is opened up for 
the low- lying antiferromagnetic fluctuations, with the weight being pushed up and concentrated at a resonancclikc 
peak at E g around Qaf = (7r, 7r), and the spin correlation length gets "truncated" at a finite scale comparable to 
the average hole-hole distance. A very unique feature in this regime is the presence of a lot of spontaneous spinon 
vortices, composites of isolated spinons locking with current vortices, which are responsible for the Nernst effect, 
residual diamagnetism, and spin Hall effect, etc. In this peculiar phase, the Cooper pair amplitude is finite because of 
the RVB pairing of the spinons and Bosc condensation of the holons. But the phase is disordered as the proliferation 
of the spinon vortices. Eventually, at a lower temperature T c , the binding of vortices-antivortices or the confinement 
of the spinons will finally lead to the superconducting phase coherence. It turns out that both the kinetic energy of 
the holons and superexchange energy of the spinons are generally benefited from this phase transition. The nodal 
quasiparticles also become coherent in the superconducting phase as the result of the phase coherence and spinon 
confinement. The non-BCS neutral S = 1 spin excitation is similar to the lower pseudogap phase with a reasonancelike 
structure at E g and the suppression of antiferromagnetic correlations below E g , as caused by the holon condensation. 

Such a top-down phase diagram demonstrates an amazing richness of the minimal phase string model, which covers 
almost all the interesting regimes of the high-T c cuprate superconductors, except for the very undcrdoped regime 
where the superconducting phase ends as well as the overdoped regime where the upper pseudogap terminates, which 
is question-marked in Fig. 31. 

At half filling, with the vanishing gauge fields, the minimal model does produce an AFLRO state with a very precise 
variational ground-state energy. But once away from half-filling, the holon condensation will force the opening up 
of a spin gap E g oc 5 J and thus the disappearance of the AFLRO. But this is unphysical at sufficiently low doping 
where the long-range antiferromagnetic correlations should remain dominant. Indeed, in this regime a topological Z 2 
vortex can become a low-lying excitation once being bound to a holon to form a localized composite object, known as 
a hole dipole. In this regime the correct low-energy phase string model is given by Eqs. (131) and (132), and a hole 
dipolc in Fig. 29 can be regarded as the realization of the dual holon confinement in the antiferromagnetic phase in 
contrast to the spinon confinement in the superconducting phase. 

In the overdoped regime where the spinon RVB pairing disappears at T = 0, the minimal phase string model should 
be also modified. A possibility is for the Z 2 vortex to be bound with a bosonic spinon such that two gauge fields 
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are effectively cancelled in Eq. (131), in favor of the holon condensation as well as the kinetic energy. In this way, 
the bosonic spinons will be effectively turned into fermions and a Fermi liquid state may be recovered if the bosonic 
holons remain condensed in the high-doping regime. So the phase string model may simply reduce to the slave-boson 
mean- field description at 5 > xrvb- 

Finally we remark on that throughout this paper, only the nearest neighboring hopping of the t-J model is con- 
sidered, which is related to the origin of the singular phase string effect. However, the phase string effect will get 
qualitatively modified in the presence of the next nearest neighbor hopping process. This is an important issue which 
has not been touched upon so far. Just like the phase string effect is very singular in the original t-J model, the 
next nearest neighbor hopping term will also become singular in the phase string formalism. It is thus expected to 
be important to interpret the detailed experimental measurements in the cuprates, e.g., the asymmetry in the holc- 
and electron-doped cuprates. 
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